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We investigate the influence of nucleon superfluidity on the neutrino emissivity of nonequilibrium 
P processes. Calculations of the reduction factors for direct and modified Urea processes with three 
types of nucleon superfluidity in npe matter are performed. The numerical results are given because 
the analytical solution is impossible. We find that the superfluid influence is closely related to the 
chemical departure from (3 equilibrium. For a small chemical departure, the superfluid reduction 
factor depends almost only on the gap and is hardly affected by the departure, while for a large 
enough departure, it rapidly enhances neutrino emissivity. The onset of "enchanced" emission has 
some corresponding thresholds that seem to be linked to the ratio of the energy gap to the chemical 
^"■r: ■ departure. 
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I. INTRODUCTION 



In superdense nuclear matter that constructs the core of a neutron star, a great amount of energy is carried away 
r-rV' by the neutrino. For the standard composition of superdense matter (neutrons, with an admixture of protons and 
electrons), the main neutrino emission mechanisms are powerful direct Urea reactions, 

<—| n— >p + e + T^, p + e^n + is e , (1) 

and much weaker modified Urea reactions, 

O 



N -> p + N + e + z/e, p + N + e^n + N + v e . (2) 



The direct Urea reactions can proceed only if the ratio of the proton number density to the total baryon number 
density n p /n exceeds a certain threshold value which allows simultaneous energy and momentum conservation [lj. 
Otherwise, direct Urea reactions are forbidden. The most powerful neutrino energy losses are produced by modified 
Urea reactions. 

Furthermore, the j3 reactions just mentioned bring the constituents into the state of chemical equilibrium, 
fJ-n = fJ-p + Me, which determines the relative concentrations of particles. However, if any external or macroscopic 
CNJ ■ phenomenon changes the density of a matter element, the /? equilibrium state will change, resulting in a departure 
\ from equilibrium quantified by the chemical imbalance Sfi — fi n — fi p — /j. e . Several authors have investigated many 
I ■ astrophysical situations in which the nonequilibrium /3 processes occur, such as gravitational collapse of neutron star 
Q\ | , pulsar spin-down , a hypothetical time variation of the gravitational constant @ , the existence of even rela- 
tively slow hydrodynamic flow [7] or of millisecond oscillations 8] in the neutron star interiors. The departure from 
equilibrium, Sp ^ 0, increases the phase space for nonequilibrium f3 process and strongly enhances corresponding 
• i-H , neutrino emissivity with respect to the chemical equilibrium value [2j . 

It is well known that below a certain critical temperature nucleons in neutron star matter will be in superfluid 
states. The appearance of energy gaps A„ and A p reduces the particle momentum space, which contributes to the 
reaction rates; that is, the "effective" widths of momentum space of the reacting particles near the Fermi surfaces 
become thinner. This suppresses the reaction rate and the neutrino emissivity. In /3-equilibrium matter, superfluid 
reductions in reaction rate and neutrino emissivity are dependent on the ratio of the gap to the temperature A/7 1 . 
For strong superfluidity, the reductions decay exponentially, while for intermediate values of A/T, the situation is 
more sophisticated. The detailed expressions are given in a reviewQ. 

The departure from equilibrium 8pL increases the neutrino emissivity, while the superfluid energy gaps behave in 
the opposite way. Indeed, the nonequilibrium f3 reaction with nucleon superfluidity is a new case that is worth 
studying. Reisenegger Q first analyzed the mechanism to produce a faster reaction owing to nonequilibrium with 
neutron and proton superfluidity. Using a crude model, he assumed that for Sp < A p + A„, /3 reactions were 
completely suppressed, while for 5p > A p + A„ the effect of superfluidity could be neglected. Villain and Haensel 
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TABLE I: Three types of superfluidity 



Type 


Superfluidity type 


A F(ti) 


T c /A(0) 


A 




1 1 


0.5669 


B 


3 P 2 (mj = 0) 


1/2 (l + 3cos 2 i9) 


0.8416 


C 


3 P 2 (\mj\ = 2) 


3/2 sin 2 ?? 


0.4926 



broke through Reisenegger's steplike modeling and calculated the precise reduction factors of the net reaction rates 
for nonequilibrium direct and modified Urea processes in the presence of various types of superfluidity by means 
of sophisticated numerical methods |1G |. However, neutrino emissivity is also an important quantity. We need to 
calculate the relevant neutrino emissivities, which, along with the net reaction rates, can be used to simulate the 
evolution of supcrfiuid neutron star cores that are off /3 equilibrium. Our primary goal here is to obtain the total 
neutrino emissivities of nonequilibrium /3 processes in superfluid circumstances. 

The paper is organized as follows. In Sec. II, we recall the nonequilibrium condition and basic features of nucleon 
superfluidity. In Sees. Ill and IV, we detail the numerical calculations of neutrino emissivities for Durca and Murca 
reactions and show our results. Section V includes a short summary of our conclusions and discussion of possible 
applications. 



II. DESCRIPTION OF CHEMICAL DERIVATION AND NUCLEON SUPERFLUIDITY 



We consider a simple npe model of neutron star matter. Each of the constituents is strongly degenerate, with Fermi 
momenta pp i and chemical potential /ii (i = n,p,e). They are in thermodynamic equilibrium but not necessarily 
in chemical equilibrium. In the absence of f3 equilibrium, there is a finite difference in the chemical potentials, 
5/i = fi n — /i p — /i e ^ 0, following the conventions of Ref. Q. 

Nucleon superfluidity occurs via Cooper pairing of particles owing to an attractive component of their interaction, 
with the appearance of a gap 6 in the particle energy spectrum near the Fermi level fllT ] . It is widely accepted that 
there are likely three types of nucleon superfluidity in npe matter, 1 So, 3 P2 (mj = 0) and 3 P 2 {\mj\ = 2), denoted 
A, B, and C, respectively (Table QJ 12]. We introduce two dimensionless quantities, £ = Su/T and v = A(T)/T, 
to describe the chemical deviation and gap amplitude, for convenience, and use natural units with h = ks = c = 1 
throughout this paper. 

To guarantee the integrity of this paper, we recall fundamental properties and characteristic quantities [l3l - [l5j . The 
onset of superfluidity is accompanied by the appearance of the energy gap 5. Near the Fermi surface {\p~Pf \ <pf), 
we have 



\/S 2 + n 2 at p < p F , e = fi + yjb 2 +rj 2 at p > p F . (3) 



Here n = v-p(p — p F ), p-p and are the Fermi momentum and Fermi velocity, respectively; fi is the chemical potential; 
and 5 2 — A 2 (T)F(-d), where A(T) is an amplitude that determines the temperature dependence of the gap and F($) 
describes the dependence of the gap on the angle i? between the quantization axis and the particle momentum. The 
quantities A and F are determined by the superfluidity type (TableU]). I n case A the gap is isotropic, and S = A(T). 
In cases B and C, the gap depends on ■&. Note that in case C the gap vanishes at the poles of the Fermi sphere at 
any temperature: Fc(0) = Fc(n) = 0. 
The gap amplitude A(T) is determined by the BCS theory which can be written as 

l^W¥f-/W (4) 



A(T) J 4tt 

where Ao = A(0), dfl is a solid angle element in the direction of particle momentum p, / = (1 + e 2 )^ 1 is the 
Fermi-Dirac distribution, A is a numeric coefficient (Table and 



sign(x)y/x 2 +y 2 , x = — , y = - (5) 

Note that in the absence of nucleon superfluidity, z = x. According to Table HI the dimensionless gap y can be written 
as: 



VA = va, Vb = v B yl + 3 cos 2 y c = v c sintf. (6) 
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In the following calculations that involve nucleon superfluidity, we assume that only one type of nucleon is sirrjerfluid, 
which is equivalent to assuming that the larger gap prevails. This approximation seems to be quite reasonable [ill ] . Our 
goal in this paper is to investigate the influence of nucleon superfluidity on the neutrino emissivity of nonequilibrium 
(3 processes, which is described by the reduction factor R x (^,Vj) — Qx(£,,Vj)/Qx(0 = lx(.^ v j)/^x(0- Qx(£,Vj) 
and Qx{£,) refer to the neutrino emissivities in two cases: one is nonequilibrium j3 process with nucleon superfluidity; 
the other is nonequilibrium /3 process without nucleon superfluidity. Furthermore, X labels the type of reactions, i 
denotes the superfluidity type (A, B, or C), and j the type of superfluid nucleon {n or p) with gap Vj. Note that for 
T > T c (the critical temperature of the superfluid nucleon), R x (^ 7 Vj) = R x (£, 0) = 1, and for T < T c , R x (£,, Vj) < 1. 

III. DIRECT URCA PROCESSES 

As shown in Ref. [l|, the direct Urea process is allowed by the momentum conservation when pp n < pp + Pf c - For 
pure npe matter where pf p = PF a , it corresponds to n p /n > 1/9. This happens if the density is several times higher 
than the standard nuclear matter density, po — 2.8 x 10 14 g cm~ 3 . 

In the absence of /3 equilibrium, the neutrino emissivities of two direct Urea processes with nucleon superfluidity 
are in the following forms: 



47T 

(270 s 



dxi f{zi) 



3 f 3 poo 

TT / dty 5(P f - P l )|M /4 | 2 TT p Fj m* / dx u xl 



dx 2 (1 - f(z 2 )) 



+ 00 



dxz (1 - f(x 3 ))6(z 2 + x 3 - zi + x v + C)i 



(7) 



Sp 



47T 



+ 00 



dx 2 f(z 2 ) 



+ ■00 



3 f 3 p oo 

IT / dCLj S(P f - P^Mf^ 2 TT PF 3 m* / dx v x\ 



dx 3 f(x 3 ) 



+ '0O 



dxi (1 - f(z 1 ))8{z 1 - z 2 - x 3 



0- 



(8) 



where j — 1, 2, 3 corresponds to n, p, e respectively, x v — p v /T = e v /T is the dimensionless energy of the neutrino, 
£, Vj, and Zj (j = 1, 2) have been defined in Sec. HH and f(x) = (1 + e x ) _1 is the Fermi-Dirac function of nucleons and 
electron, pp s is the corresponding Fermi momentum and m* is the effective particle mass. Furthermore, dVLj is the 
solid angle element in the direction of the particle the momentum pj, and the 5 functions describe momentum and 
energy conservations of the particles in the initial and final states. Finally, |-M/i| 2 is the squared reaction amplitude, 
and for nonrelativistic nucleons, it is independent of particle momenta and can be taken out of the integral. In this 
paper, we focus on the total neutrino emissivity Qd = Qn + Qp°^- 

Let us start with the total neutrino emissivity Qq of the direct Urea processes without nucleon superfluidity under 
j3 equilibrium. Under the condition of j3 equilibrium and without nucleon superfluidity, we set £ = and replace Zj 

with Xj . The direct and inverse reactions have the same neutrino emissivity. Then Qq D ^ can be written as (for details 
see Ref. @): 



(2tt) 8 

3 



T 6 A D I Q 



A n = Air 



AD) 



3 

IL 

3=1 



3 = 1 J 



<KP/-Pi), 



n 

3=1 



+ OC: 



dx 3 fj 



/ j -^3 
J =1 



(9) 
(10) 

(11) 



Here the integrals Af> and Iq D ^ are standard ([l 



An = 



32tt 2 



VF n PF v PF^ 



(D) _ 



4577r b 
5040 



(12) 
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A. Nonequilibrium without superfluidity 



In this section, we recall the neutrino emissivity of nonequilibrium Durca processes without nuclcon superfluidity. 
When /3 processes are off chemical equilibrium, £ =^ 0, the direct and inverse reactions of Urea processes have different 
neutrino emissivities. The expressions for the neutrino emissivities of the two reactions are similar to Eqs. (J7J and 
([8]) , but the variable Zj is replaced with the corresponding Xj . The total neutrino emissivity can be presented in the 
forms @, 1 



with 



According to Ref. Q, 



We obtain 



Qd(0 = j^T 6 A D I D (0\M fi \ 2 f[p Fj m*, 



/ D (0 = / dx v x% [J{x v - + J{x v + 0] 
Jo 



J{x v -0 = 

J{X V +0 = 



:>> 



n 

3 

n 



dxj fj 



3 = 1 ' 
3 /»+oo 



dxj fj 



S(xi + x 2 + x 3 - x v + £), 



*i>(0 = i 



1071£ 2 315f 4 21£ 6 



457tt 2 457tt 4 457tt 6 ' 



^(0 



45 77r b 
2520 



1071£ 2 315£ 4 21£ 6 \ 
1 + 457tt 2 + 457tt 4 + 457?r 6 J 



(13) 

(14) 
(15) 

(16) 

(17) 
(18) 

(19) 



B. Nonequilibrium with superfluidity 

Consider suppression of nonequilibrium Direct Urea processes by proton or neutron superfluidity. The superfluidity 
affects the dispersion relation of nucleons under the integrals in Eqs. (fT5|) and (fT6|) . in accordance with Eq. (|3|). The 
neutrino emissivity can be written as 



(20) 



Vj) is the superfluid reduction factor. The total neutrino emissivity Qd(£,Vj) = Qn (£, ^j) + Qp D \£,,Vj) is in 
the following form: 



Qd{Lv 3 ) = j^T 6 A D I i D (Z,v j )\M fi \ 2 np Fj m* j , 



(21) 



where 



4tt 



A 



D 



3 

n 



S(P f - Pi) / 



+ 00 



dxidx 2 dxz [A(l - /2XI - /3)<5(^i - zi - x 3 - x v + f) + ^2/3(1 - fi)S{zi + x v - z 2 - x 3 + 0] > ( 22 ) 



and dfij is the element of the solid angle in the direction of the momentum of superfluid nucleon j. 
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Here, we consider the case where only one type of nucleon is superfluid, then the z variable for the nonsuperfhiid 
nucleon has to be replaced with the corresponding x variable, while Vj is "included" in the z variable for the superfluid 
nucleon. Thus, the reduction factor of nucleon superfluidity for a nonequilibrium Durca process can be presented by 



MO 



i 



MO 







r/2 

BmWdtiHhi&Vj), (23) 



where 



poo r+oo 

H L>(€> v j)= / dx u xl I dx 1 dx 2 dx 3 [/i(l - / 2 )(1 - fz)5{zi ~ z 2 - x 3 - x v + £) 

JO J -oo 

+/ 2 / 3 (l - fi)S(zi +x„~z 2 -x 3 + 0]. (24) 

In Eq. ()24|) . integration over the electron variable can be done using the formula for Fermi integrals, similar to the 
treatment in Rcf. 10], 



dxf(x)f(y -x) = -JL_=G(y) (25) 



- DG 



which enables us to integrate over the x variable for the nonsuperfluid nucleon, giving 

oo +oo 

#£>(£>%) = / dx » x l I dx j [f{z j )G{x v - f - Zj) + f(-Zj)G(x v -£ + Zj) 



+f( Zj )G(x u + £ - Zj) + f{-z 3 )G(x u + £ + z )]. (26) 

We assume only one type of nucleon superfluidity and numerically evaluate R l D (£,,Vj) for the wide range of interest 
with algorithmic methods (refer to the appendix of [13] )■ lu Fig- Q]we present the reduction factor for Durca reaction 
with isotropic superfluidity of protons or neutrons. The variable is dimensionless chemical departure \og 10 (Sfj,/T). 
The curves correspond to different values of dimensionless gap amplitude log 10 (gap/T). From this figure we can get 
that there are at least two features: (I) for a given curve, the change of the chemical departure value almost does 
not affect the reduction factor until it increases at a value; (II) for huge gaps log 10 (gap/T) > 1, as soon as Sfi > gap, 
and until the reduction factor becomes close to 1, the curves are almost parallel. This means that for small chemical 
departure, the reduction factor can be written as the function of the variable gap only, while for huge gaps it can be 
written as R ~ f(gap)Q(gap/5p,). These are similar to the net reaction rates presented by Villain and Haensel [lfj| . 

Additionally, the behavior of the superfluid suppression on the neutrino emissivity here is different from that 
proposed by Reisenegger He suggested that for 5fi < A p + A n the beta reactions were completely suppressed 
(R = 0), while for Sfi > A p + A„ the effect of superfluidity could be neglected (i? = 1). However, as can be seen 
from this figure, the reduction factor R doesn't equal to or 1. In other words, we find that in a proper range 
of the value Sfx/T, the neutrino emissivity will be suppressed by the nucleon superfluidity, but it is not completely 
suppressed. If taking the onset of rising R as the threshold point, we find that in the region of \og 1Q (gap/T) < 1.0, 
\og w (5 [ith / T) = constant (approximately 0.78 in Fig. [5]), while for \og w (gap/T) > 1.0, 6/j,th ~ gap. Here fith 
is the value of p, at the threshold point. Eventually, when and only when the value of gap/T approaches infinity, 
gap/T — s> oo, our results are fully same with that in Reisenegger jj]. These are easily illustrated by plotting K as 
function of Sp,/T in Fig. (5J where K is the partial derivative of R with respect to \og 1Q (Sp,/T). It is clearly that the 
curves in this figure can confirm our above discussions. 

The above behavior of superfluid suppression on the neutrino emissivity can also be seen from all the other com- 
parable figures in the following. 

Figs. [3] and S] show the dependence of the reduction factor of Durca reaction on dimensionless chemical departure 
and dimensionless gap amplitude with anisotropic superfluidity of B and C, respectively. Compared with Fig. [TJ we 
can see that the reduction factors can either increase (type C) or decrease (type B) in respect to that of type A, which 
means that the maximal value of the gap on the Fermi surface is important since the value for case C is larger than 
for case A, while for case B is smaller. 



-4-3-2-101234 

log 10 ( Su/T) 



FIG. 1: Dependence of the reduction factor on the chemical imbalance for various gap's amplitude. The numbers next to the 
curves give value of \og 10 (gap/T). Direct reaction with isotropic superfluidity of protons or neutrons. 




FIG. 2: Dependence of K (K = dR/d(log 10 (5jj,/T))) on the chemical imbalance for various gap's amplitude. The numbers 
next to the curves give value of log 10 (<7ap/T). Direct reaction with isotropic superfluidity of protons or neutrons. 



IV. MODIFIED URCA PROCESSES 



As mentioned in the previous section, the direct Urea process in npe matter is allowed at densities several times the 
nuclear density. Then in smaller densities the modified Urea process becomes the main neutrino reaction. Reaction 
([2]) differs from reaction ([T]) by the presence of an additional nucleon to ensure momentum conservation. The modified 
Urea process will be labeled by upperscripts (MN), where N = n indicates the neutron branch of reaction @ and 
N = p indicates the proton branch. Similar to above section, we give the expressions of the neutrino emissivities for 
the direct and inverse Murca processes out of equilibrium in superfluid nuclear matter, 

5 /.OO 

5^ f -Pi)\M fi \ 2 J[p Fj m* / dx v 4 

3 = 1 J ° 

/+oo 
dzida^dxsda^dzs f(z 1 )f(z 2 ) (1 - f(z 3 )) (1 - /(^4))(1 - f(x 5 ))5(z3 + z 4 + x 5 - Z\ - z 6 + x v + £), (27) 
-oo 
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FIG. 3: Same as in Fig. {TJ but with anisotropic superfluidity type B of the neutrons. 
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FIG. 4: Same as in Fig. fTJ) but with anisotropic superfluidity type C of the neutrons. 
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On 5 />OG 

1 = 1 J 4=1 J ° 



dccida^dzsda^da^ (1 - f(z\)){l - f(z 2 ))f(z 3 ) f(z4)f(x 5 )S(z 3 + z± + x 5 - zi - z 2 - x„ - £). (28) 



where j — 1, 2, 3, 4 correspond to n, Ni,p, Nf (iVj and Nt refer to the initial and finial spectator nucleons), while j = 5 
corresponds to the electron. The remaining parameters are the same as in Durca case. The total neutrino emissivity 
of one branch of Murca processes is Q^ MN ^ = Qn + Qp MN ^ . 

If none of the nucleons is superfluid, we have no necessary to make distinction between the two branches in the 
phase-space integral and Zj — Xj respectively. Under the condition of chemical equilibrium, £ = 0, the neutrino 
emissivity of direct and inverse Murca processes without nucleon superfluidity are the same. The total neutrino 
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emissivity of one branch of Murca processes have been calculated by [TtJ , 



n (MN) _ 1 T 8 a t(MN) , m ,2TT„ , 

~ f27r) 14 Amn1 o \ M fi\ llPFj 1 

1 ' 3 = 1 



Amn = 47r 



AMN) _ 
'0 — 



5 

n y d ^ 



3=1 



d^iy X,, 



n 

3=1 



3=1 



(29) 
(30) 
(31) 



Here fj = (1 + e Xj ) 1 is the Fermi-Dirac distribution function. The quantity |Mfi| 2 appear to be constant and can 
be taken out the integral. The detailed calculations refer to @. we have 



Amti — 



27r(47r) 

Pk 



2(2tt) 5 (mat) H513^ 8 
Am p = ~ ^3 ~ (Pfc + 3^ " PfJi h ~ - 



PfJ>%Pf c 



60480 



(32) 



A. Non-equilibrium without superfluidity 

The neutrino emissivity of non-equilibrium modified Urea process without nucleon superfluidity has been studied 
by some authors (@, Q)- When the matter is out of beta equilibrium, <5/i ^ 0, it opens additional volume in the 
phase-space which results in an increase of the neutrino emissivity. It can be written as: 



Qmn{£) 



1 5 
2^)T4 r 8 AmnIm n (0\Mfi\ 2 I] P*i K > 



with 
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J o 
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Then the effect of the equilibrium departure on the neutrino emissivity of Murca process can be presented as: 

Qmn(0 _ Imn(Q 
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SO 



21, 



(MN) ' 



(I 



According to 



(0 = Q ( o MN) F M ({) 

22020£ 2 5670£ 4 



F M (0 = 1 + 



420<f 



9£ 8 



115137T 2 115137T 4 115137T 6 115137T 8 



Thus 



Imn (0 



420£ 6 9£ 8 \ 

115137T 2 ' 115137T 4 ' 115137T 6 + 115137T 8 J 



JIVCS;,' 4 ( i 22020£ 2 , 5670£ 4 
60480 ' ' 



(33) 

(34) 
(35) 

(36) 

(37) 

(38) 
(39) 

(40) 
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B. Non-equilibrium with superfluidity 



The analysis of the modified Urea process is quit similar with the direct Urea process. The appearance of nucleon 
superfluidity reduces the available phase-space which suppresses the neutrino emissivity, Zj ^ Xj. We label R l MN (£,i Vj) 
as the superfluid reduction factor. We can write 



(41) 



Notice that for T > T c (the critical temperature of superfluid nucleon) R\ IN (^, Vj) = R 1 mn(£,, 0) = 1, and for T < T c 
& MN (£,Vj) < 1. The total neutrino emissivity Q( MN \£,Vj) = Q ( n MN) (tvj) + Q P MN) {Z,Vj) is in the following form: 



Qmn(Z,Vj) = 

where Amn is the same with Eq. (|30[) and 

47T 



2(2tt) 
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PFjl 



(42) 
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MN 



n / ^ s(p f - pa / a Xv xi n / **i i 
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5 5 

[6{J2 Zj -Xv + Z) + *j -x v - 01- 

3=1 J'=l 

(43) 



Thus, the reduction factor of nucleon superfluidity for non-equilibrium Murca process can be presented by 



t-MN 



(0 ' 



(44) 



In the following, we numerically analyze some special cases that may occur in dense matter and give the diagrams 
to illustrate the relationship of superfluid reduction factor with the equilibrium departure and the superfluid gap 
amplitude. 

Case 1. the neutron branch with isotropic superfluidity of protons 

Since the singlet-state pairing proton gap is isotropic, the angular and energy integrations in Eq. (|43|) are separated. 
The angular integrals are the same as in the non-superfluid case, and we get 



hln (£) 



1 



'Mi 



with 



1 f°° 
Hi(x) = - dss' 
6 Jo 



exp(s — x) — 1 



dx p f{z p ){H l {z p + £) + Ht(z p - £)), 



[(s - xf + 4tt 2 



(45) 



(46) 



This result is also valid for the case when the neutron superfluidity is of singlet type in the proton branch of Murca 
processes. 

Case 2. the proton branch with any superfluidity of neutrons 

Let the protons be normal, but the neutrons be superfluid due to nn Cooper pairing. The integration over solid 
angles dClj of all particles but the neutron in Eq. (l43l) is the same as for non-superfluid case. We get 



ft nX (f „ ^ - I Mp{£,,Vn) 



1 



Iup{£) Im p (Q Jo 



dXnftZnXH^Zn + f) + Hi(z n - £)), 



(47) 



where <; = costf, ?? is the angle between the quantization axis and the momentum of the superfluid neutron and Hi(x) 
is the same as in Eq. (|4"5j) . 

Case 3. the neutron branch with isotropic superfluidity of neutrons 

Now consider the neutron branch with isotropic superfluidity of neutrons. Let j = 1,2,3 in Eq. (|43[) refer to 
neutron. Integrating over X4 and x 5 , we get 



'Mi 



wo 



1 3 r + co 

wiyny^ ■ 



dxjf{zj){Hi{zi + z 2 + z 3 +£) + H 2 (z 1 +z 2 + z 3 - £)), 



(48) 
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FIG. 5: Dependence of the reduction factor on the chemical imbalance for various gap amplitude. The numbers next to the 
curves give value of \og 10 (gap/T). Modified neutron branch with isotropic superfluidity of protons. 
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FIG. 6: Same as in Fig. ([5]) but for Modified neutron branch with isotropic superfluidity of neutrons. 



poo 

H 2 (x) = / dss 3 S ~ X (49) 
J exp(s - x) - 1 

Employing the numerical calculations similar with that done in we give the results for Case 1 and Case 3 
in Figs. [5] and [6] respectively. We need not give the similar graph for Case 2. The result of Case 2 with isotropic 
superfluidity of neutrons are same with Case 1. Comparing Figs. [S] and [51 we can see that the impact with superfluid 
spectator nucleons is much larger than that with superfluid non-spectator nucleons. Based on the similar reasons 
with that presented in [ifj| . the cases with triplet-state nucleon superfluidity have been neglected. 



V. DISCUSSION AND CONCLUSION 

We apply the methods of Villain and Haensel in paper [lO] to calculate the effects of nucleon superfluidities on the 
neutrino emissivities of direct and modified Urea processes out of equilibrium, while they calculated the corresponding 
net reaction rates. Although the numerical methods are similar, there are two differences in calculating neutrino 
emissivities and net reaction rates which can be seen easily by camparing Eq. (|43[) and Eq. (25) of paper [Toj . for 
example. Firstly, the integral / dx u x z v takes the place of / Ax v x 2 v . Secondly, the two delta functions add together 
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in Eq. (l43t while subtract each other in Eq. (25) of previous paper. So the calculations of the neutrino emissivities 
are necessary. We here give the graphical overview of those results for 

- Durca reaction with superfluidity of each type; 

- Murca reaction with isotropic superfluidity of the non-spectator nucleon (e.g. proton for the neutron branch); 

- Murca reaction with isotropic superfluidity of the spectator nucleon (e.g. neutron for the neutron branch). 

As shown in calculations, the appearance of nucleon superfluid reduces the neutrino emissivity by a factor that 
depends on the value of A/T and of 8fi/T, the types of the beta process (direct or modified Urea) and of superfluidity. 
When the quantity 5/i/T is relatively small, the superfluid reduction factors don't change with the value of Sfi/T and 
become the function of only one variable A/T. It should be similar with that under the condition of beta equilibrium. 
However, the violation of superfluid suppression on neutrino emissivity occurs when the chemical departure reaches 
a threshold value. The threshold tends to be the result presented by Reisenegger (Q) only when the energy gap of 
superfluidity is large enough. 

It's well-known that neutrino emissivity plays a key role on the cooling of neutron stars ([H1[I(|). we recognize that 
although the superfluid gap suppresses dramatically the neutrino emissivity in dense matter, it is possible that the 
chemical imbalance of beta reaction large enough will break the suppression. Meanwhile, the imbalance contributes 
chemical heat to neutron star core. Along with the results of reaction rates calculated by Villain and Haensel ((Io|). 
our results are required for numerical simulation of neutron star cooling. We have the primary investigations on this 
problem. It will be reported in our another paper. 
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